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A general theory for the core-level electron excitation of anisotropic systems using angular inte- 
grated electron energy-loss spectroscopy has been derived. We show that it is possible to define a 
magic angle condition at which the specimen orientation has no effect on the electron energy-loss 
spectra. We have not only resolved the existing discrepancy between different studies of the magic 
angle condition, but also extended its applicability to all anisotropic systems. We have demonstrated 
that magic angle electron energy loss spectroscopy is equivalent to the orientation averaged EELS, 
although the specimen remains stationary. Our analysis provides the theoretical framework for the 
comparison between theoretical calculation and experimental measurement of core-level electron 
excitation spectra in anisotropic systems. In addition to MAEELS, we have also discovered a magic 
orientation condition which will also give rise to orientationally averaged spectra. It's relation with 
the magic angle X-ray absorption spectroscopy and magic angle spinning nuclear magnetic resonance 
is demonstrated. 

I. INTRODUCTION 

Anisotropic systems differ from isotropic ones in that their response to an applied force or field depends not only 
on the magnitude but also the orientation of these influences^. For electron energy loss spectroscopy (EELS) or 
X-ray absorption spectroscopy (XAS), a good example is the carbon Is core electron excitation in graphite2iiMi£i£i 
Excitation into the unoccupied states of 7r-symmetry is only allowed if the applied field is along the direction normal 
to the graphite sheet (defined to be the local z-axis). Excitation into the unoccupied state of cr-symmetry can only 
occur if the applied field lies in the plane of the graphite sheet (defined to be the local x-y plane). As a consequence, 
the intensities of these two excitations depend on the specimen orientation in general. Graphite is an uniaxial system 
which is the simplest example of anisotropic systems. Anisotropic response is a gift to the experimentalists as it offers 
further insight into the electronic excitation process^i^ or it can be used to determine the orientation of molecules&iS 
or internal magnetic fieldii. In this paper, we focus on the complexity it brings to the EELS measurement which now 
also depends on the precise orientation of the specimen and ways to overcome it. 

Concerning electronic excitation, many important systems show anisotropy such as familiar layered materials 
like graphite and BNiS, non-central symmetric semiconducting compounds such as GaNA£, superconductors such 
as YBaCuGi2*i£*i£ii£ and MgB2 1819-20 ' 21 . Some nanostructures such as nanotubes and nanoonions can be considered 
as roll-up versions of layered materials, so the local anisotropy also changes with position inside the nanostructured 
material 2 ^. In addition, shape or local field effects will turn an isotropic transition into an anisotropic one^i. This 
presents difficulty for quantitative spectral analysis, particularly from localized areas using EELS, a powerful method 
for electronic structural characterization, particularly at nanoscale. For example, real changes in the electronic struc- 
ture may not be easily distinguished from variations due to mere specimen rotation or probe displacement along a 
curved basal plane. So there is a need to find an experimental approach in which specimen orientation is no longer 
effective on the spectra. 

A related problem also exists in other forms of spectroscopy. For instance, in solid state nuclear magnetic resonance 
(NMR), the signal depends strongly on the specimen orientation and more useful information is obtained after the 
discovery of the magic angle spinning nuclear magnetic resonance (MAS-NMR) methocUi* 2 ^^. Another example in 
XAS is the so-called 'magic angle' condition used to eliminate the orientation effect for anisotropic systems^. In 
EELS, the applied field responsible for electronic excitation is parallel to the direction of the momentum transfer 
vector q of the incident electron in the inelastic scattering 8 . As shown in fig^ by virtue of electron scattering, the 
direction of q is a function of the electron scattering angle 9, being parallel to the incident electron beam at the zero 
scattering angle and changing towards a direction perpendicular to the beam with the increasing scattering angle. 

In the angular integrated energy-loss spectroscopy with a centered circular aperture, a standard EELS collec- 
tion condition particularly for high spatial resolution and high signal to noise ratio, electronic excitations of dif- 
ferent momentum transfer q are recorded. Excitations induced by the applied fields over a range of directions are 
summed22i2L2Si2S. Menon and YuanSi showed both experimentally and theoretically there also exists a magic angle 
(MA) at which the EELS signal is independent of specimen orientation. Their theory predicts the magic angle for 
a parallel beam illumination (MAm) to be 46e- There 9e is the characteristic angle which has a approximate form 
6>£=2E/Eq (the more accurate definition is given by Ritchie and HowieSi) and with E being the energy-loss and Eg 
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the incident energy of the fast electrons. The derivation is for uniaxial, graphite-like, systems. 

A number of other studies have also analyzed this problem again, mostly for uniaxial anisotropic systems. However, 
controversy as to the precise magic angle condition persists. For example, Zhu et alm^ used an approximate method 
to arrive at a value of MAm of about 1.89e- Paxton et al— used a different approach to arrive at MA||=1.36 
9e- Compounding this state of confusion in magic angle prediction is the report of Daniels et al— who found 
experimentally that MA||=2#£; and also provided a theoretical justification for their result. Souche et aliSS, presented 
a more through theoretical analysis of the anisotropic electron energy-loss spectroscopy for graphite-like materials, 
taking into account the convergence beam effect. Their result for the parallel beam illumination is in agreement with 
the prediction of Menon and YuanSi. In addition, apart from attempts by Paxton, all theoretical analysis has been 
confined to MAEELS in uniaxial systems, the applicability of the magic angle concept in more complex anisotropic 
systems is an unexplored area. 

In order to make use of MAEELS, it is vital to understand the discrepancies between the various theories, par- 
ticularly given that all models make the same basic assumptions, i.e. the electron beam scattering being kinetic, 
non-relativistic, and obeying the dipole approximation. In the process, we have demonstrated that the magic angle 
is a general effect which can occur in all sorts of anisotropic systems, be they crystalline, amorphous, single phase 
or nanostructured and applies to anisotropic transition whether it is intrinsic or extrinsic (for example, shape-effect 
induced). 

The paper is organized as follows. In section II, we present a through analysis of electron energy-loss spectroscopy in 
anisotropic systems in order to establish a solid foundation for MAEELS in the core-electron excitation region and we 
derive a general definition for magic angle conditions. In section III, we show that the spectral information obtained 
at the magic angle is equivalent to orientationally averaged spectroscopy. From this geometrical interpretation, 
we can understand the different approaches used in MAEELS analysis reported in the literature and explain the 
reasons for discrepancies in the prediction of magic angle conditions. In section IV, we have provided more explicit 
MAEELS expressions for applications in specific material systems. In particular, we concentrated on more commonly 
encountered uniaxial and biaxial systems and establish the connection between MAEELS and other forms of magic 
angle spectroscopies such as XAS and MAS-NMR. 

II. GENERALIZED MAEELS THEORY 

There are two ways to study the electronic excitation in anisotropic systems: the macroscopic approach is to use 
the dielectric function and the microscopic approach to calculate directly the quantum mechanical transition matrix 
element. Both approaches have been reported in the literaturei^22i2§i2L2ii22i2ii and physically they should yield the 
same result. To facilitate comparison of these reports, we will make our derivation using both approaches. It is 
known that the quantum mechanical approach can reveal the microscopic physics involved, but because of the need 
for accurate wave functions, it does not always yield useful experimental results. On the other hand, the dielectric 
approach is a phenomenological description of materials response that can be measured accurately, even though the 
microscopic origin of the electronic transitions may be obscured— 

A. Dielectric Formalism 

Here we start with the dielectric approach where the calculation is relatively straightforward, because the informa- 
tion required is not the detailed excitation form of but the overall effect in terms of the response to a perturbation 
deserved by the relation^: 

n, Y, " "'v (!) 
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where D, the electric displacement, is related to the electric field E by the well-known e 4 - 7 the dielectric function of 
the material system, e 4 -? is a 'metric tensor—, defined in terms of a reference frame where the orthogonal principle 
axes are aligned with the major symmetry directions of the physical system. For convenience, we have defined this 
reference frame as the sample frame (x,y,z). 

The imaginary parts of (-1/e) is known as the energy-loss function&SiSLSi, providing a complete description of the 
response of the medium through which the fast electron is travelling. The double differential cross section used to 
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estimate the intensity of EELS in an anisotropic system can be expressed as&21: 

da 2 (q) _ Am 



dEdtt na a h 2 V)gj£ ij 'gj 



Im{-——) (2) 



where m e the mass of the fast electron and q is the momentum transfer vector of the fast electrons in the inelastic 
scattering process (see fig^), n the number of atoms per unit volume of the material, a® the Bohr atomic radius, h 
the Plank constant and q; the projection of q in the sample frame. Note the components of the dielectric function 
are assumed not to be a function of q, and this assumption corresponds to the dipole approximation in the quantum 
mechanical analysis of single electron transitions 3 ^. 

We will restrict our discussion to core electron excitations in which we use the approximation&SSiSl for ei=Re(e) 
and £2=Im(e), e\ ks 1 and £2 ~ 0. Eq.(J2| can then be simplified to 

da 2 (q) _ 4m e \ ^ q^£ 2 tJ g 3 

dEdn ~ na h 2 ^ q 4 U 

1,3 

From now on, we will be only interested in the EELS which is obtained by integrating Eq. over the angular range 
determined by the collection condition, i.e. the convergence semi-angle ao for a convergent beam and the collection 
semi-angle (3q for the centered collection aperture used. This gives the partial angular integrated cross section 

%{a^M = ±^J2([ d^W* = ^^Y, W ^o,f3o,d)e 2 v (4) 
dhj nanh 1 ^— ' / o 4 nanh^kn ^-^ 

hJ »,J 

where O denotes the orientation of the sample with respect to the electron beam direction. Here the weighting factor 
depends both on the specimen orientation and the experimental condition used, i.e. ao and 0o. In order to 
separate out these two effects, we have transformed the representation of q from the (x,y,z) orthogonal coordinate of 
the sample frame to the (X,Y,Z) orthogonal coordinate of the laboratory frame, with the Z axis defined to be optical 
axis of the electron beam. We denote the components of q in the (X,Y,Z) frame as q'i(fig|2J)- Two representations of 
q are related to each other through the rotational transformation matrix R as: 



(5) 



For transformation between two orthogonal coordinate systems, the rotational matrix element is defined to be 
the direction cosine between the basis vector e' in the (X,Y,Z) frame and the basis vector e in the (x, y, z) fram e 1 ! 35 : 

Rij = ej • e'j (6) 

Substituting Eq.JSJ into Eq.Q, we get the definition for the weighting factor as: 

.2 



(7) 



In this way we have successfully separated out the orientation factors, in the form of the product of matrix elements, 
from the integral within the bracket which is solely determined by the experimental set-up. By inspection, we can see 
that the integration over the full azimuthal angle of vector q of the integrand with the cross-indices vanishes because 
of the rotational symmetry This means that the integral has the simplified forms as follows: 

f q i 2ir 

, , J dn ~[ = ^rf||(ao,A>) (m = n = 3) 

J d ^ = U/^4 = ^(00, ft)(m = n = lor2) (8) 
(m 7^ n) 

where we have introduced notation qy and q_L to denote the components of q that are parallel and perpendicular to 
the incident beam direction respectively (figQ and k is the magnitude of the wave- vector for the fast electron beam. 
2ir 

The factor — j is used to eliminate the dimension of the reduced integral variable £y and 
k 
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Putting the integral expression in Eq.JSJ back into Eq.(J7J) the weighting factor can be written as: 

Wij(a o ,0o, O) = ^R 3i R Sj + ti±(RuRi 3 + R2iR 2j ) (9) 

If we rearrange the product of the matrix element by applying the orthogonal property of the transformation 
matrbii as ^2 R m iRmj = ^ij, and through explicit calculation using Eq.©, we can obtain a more revealing definition 

TCI 

for the weighting factor 

Wij(a ,(3 0> O) = £x<% + - £x)cosXiCosxj (10) 

where the Xi is the angle between the optical axis (the Z-direction) of the laboratory frame and the ith basis vector in 
the sample frame. It is clear that the second term in Eci. (|10|) gives the information about orientation of the sample, 
so the magic angle condition is defined by the relation 

e l„MA oMA\ c r„MA oMA\ /i i \ 

Note that our derivation has not exploited any special symmetry properties of the dielectric function, for exam- 
ple, some specific symmetry property brought by a crystal structure. Hence the magic angle condition is valid for 
all anisotropic systems, i.e., it not only applies to single crystals, but also to amorphous materials, powders, or 
nanostructures as long as an effective dielectric function tensor can be defined and measured. 



B. Quantum Mechanical Theory 



The dielectric function approach is very useful in treating practical problems. However, to understand the micro- 
scopic origin of the electronic transition responsible, it is better to work explicitly in terms of a quantum mechanical 
theory. With the wide spread use of ab-initio quantum mechanical calculation methods, proper treatment of the 
electronic excitation will become be routine^2L2£. It is vital to know how to relate them to measurements where 
specimen-orientation is an additional variable. 

IN quantum mechanics the inelastic scattering of high energy electrons can be described adequately by the first 
Born approximation assail: 

where ao is Bohr atomic radius, and vector r is the coordinate of the electrons in the sample, the initial and final 
states of which are represented by (i| and (/| , respectively. Since we have not considered the screening of the fast 
electron Coulomb potential by other electrons inside the material, this expression is only applicable to core electron 
excitations. Because of the inverse q-dependence, electron scattering is concentrated at small angles. We can then 
expand the matrix element in terms of q and only retain the first non-zero term which is the dipole approximation 8 , 
to obtain 

After projecting q in the sample frame (x,y,z) mentioned above, one can see the connection between Eq. (|13fl and 
Eq.@ by identifying Im{e l °) cx (xi)(xj)* . The rest of the derivation can follow the procedure used in the dielectric 
formalism. Thus, we should arrive at the same conclusion as Eq.(JZJ), so the magic angle condition should be the same 
as Eq.jnj, i.e. £|| = £±. 



C. The solution of the Magic angle condition 



The magic angle condition refers to the convergence and collection semi-angles, ao and (3q respectively, which define 
the experimental set-up where Ea. 1)11(1 is satisfied. We recall that the fast electron has the simple energy- momentum 
fi 2 k 2 

relation Eq = - so the energy-loss process must satisfy the following energy and momentum relations 



2m 



E = E-E = tL (14) 

Ira 
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k - k f (15) 



q 2 = k 2 + k) - 2k k f cos 9 (16) 

The simplest test case for magic angle conditions is for an experimental set-up involving parallel beam illumination. 
In this case, the scattering angle involved (9) is just the function of the semi-angle (/?), defined to be the angle between 
the wave vector of the scattered electrons and the electron optical axis. For an axially placed circular detector, the 
maximum and minimum values of the momentum transfer are given by: 



fam =k -kf = k 9 E (17a) 
V 



q-max = k 2 + kj~ 2k k f cos/3 (17b) 



Using above expressions for q to calculate the integral, following Paxton et ali^i, one obtains the result for the reduced 
integrals defined in Eq. (JHJ as: 



£|| =A = -JT~ )+oTrlii- + - — — (g maa; -g mm ) (18a) 



-i-mtn Vmax 
7 max 



e± = ^ = kln^-£||) (18b) 

This complex solution for the parallel beam illumination can be simplified because we are interested in the small- 
angle region where dipole approximation holds, so we can use the small angle approximation for qx (~ k$9) and qy 
(~ Qmin = fco^e) as to obtain a more simplified form for the reduced integrals 

fii = A « — ^ (19a) 

" 2(/3 2 + l) 

< - B - A ~ 1 [H 1 + P 2 )-Jl-} (19b) 
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where we have used the reduced collection angle (3$ — Pq/Oe- We can solve the magic angle condition using the magic 
angle relation £|| = £_l i.e. B = 3A. Within the small angle approximation, this is satisfied for (3q IA = 3.970b, or 
A9e approximately as shown initially for a uniaxial systemSi. In fig|21 the magic angle solution for the parallel beam 
illumination set-up is plotted as a function of 9e- The small-angle solution is found to be acceptable for the normal 
energy-loss and collection condition as the deviation from the more exact solution obtained from Ea. fl8p and l|19|) is 
less than 5 percent. 



D. Convergence angle effect 



In many cases, explicitly for a focused probe system or implicitly because of the need to increase the illumination 
level at the sample using a slightly convergent beam, the magic angle solution for the parallel illumination condition 
becomes inapplicable. To take into account the convergence effect, we need to reexamine the momentum conservation 
relation shown in Ea. H16|) . This vector relation can be decomposed according to the vector components parallel and 
perpendicular to the electron optical axis. In the small-angle approximation (see figEJ, the parallel version of this 
relation is a scala equation: 

?„ = -q z = k 9 E (20) 
and perpendicular components of q are defined as follows: 

{qx = ko (a cos 4> — (3 cos tp) 
qy = ko (a sin (f) — (3 sin tp) (21) 
ql = q X +<lY = k 2 ° 2 

where <f> and tp are azimuth angles for wave-vectors of the incident and scattered electrons respectively and 9 is the 
scattering angle between the direction of the incident electron and that of the scattered electron. The latter is related 
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to a ((3), the angle between the incident (scattered) electron and the electron beam axis, by spherical trigonometry: (see 
appendix) 

9 2 {a, 0,<f> - ip) = a 2 + 2 - 2a/?cos(0 - tp) (22) 
Then, the reduced integral defined in Eq.(JHJ has the following form: 

f || = A' (23a) 

R' — A' 

U = (23b) 



where 



A' = — / ada / d(j) / 0d0 / d^— 2 — E a — (24a) 

2* 2 a 2 Q J Jo Jo Jo r [9 2 E + 6 2 (a,0,<t>-ip)] 2 

1 f<»0 /"Sir /-/3 <.27r ^ 

B' = — / ada / d0 / /3d/3 / d V - 7 — — — (24b) 

WuUo Jo Jo Jo &b + ° ^ M ~ <f) 



The integral can be evaluated directly as an explicit function of the convergence and collection angles as: 

A'(ao,(3o,0 E ) = -X-iD-C) 
Aa 2 

B'(a , O , M = ^{ Wl ^ - + °~ ^ + 2a 2 In D + C ^ 2 ^ ] {D C )| 

where 



C= v /-4ag/3 2 + (l + ai + /3 2 ) 2 

5.2 , 52 



D = df5 + ^ + 1 

00 = 00 1 ®E 

Thus the solution of Ea. (|ll|) is equivalent to B' = 3A'. Ea. (|23|l is equivalent to Ea. (|19|) when the convergence angle 
ao approaches zero. The numerical result is plotted as a contour in fig|S] Our result agrees with the result given by 
Souce et al 2 ^ in their determination for an uniaxial system through more tedious integration. It is worth pointing 
out that an earlier prediction by Menon and Yuan 2 ^ for the magic angle condition in the convergent beam case is 
incorrect because it did not perform the actual azimuthal angular integration for both the incident beam and the 
scattered beam. 

A striking feature of the solution shown above is the symmetry, that is, the interchangeability between the beam 
convergence range ao of the incident electron and angular range 0q of the collection of scattered electrons. This is 
already evident in Ea. (|22|l . This can be traced to the small-angle approximation. Fig|U shows that the surface of the 
Edward sphere described by the fast electrons becomes a plane of constant energy in the small angle approximation. 
It has been argued that dependence on a and 0o are thus not symmetrical 2 ^, because the interchange of a and 
produces another momentum transfer vector which has the same parallel vector component but with a perpendicular 
component which is pointing in the opposite direction. Thus in the general case, the interchangeability of the 
convergence and the collection angles may not be valid . But in angular integrated spectroscopy, and in the small 
angle approximation, cross-section depends only on the modular square of q(see Eq.JSJ), hence the interchangeability 
holds. 



III. DISCUSSION OF THE THEORETICAL MODEL 



Our theoretical analysis not only gives us a definition of the magic angle conditions, valid for an arbitrary anisotropic 
system, but it can also give a general expression for electron energy loss spectroscopy in anisotropic systems. This 
allows us to investigate in more detail the physical meaning of the magic angle conditions as well as understanding 
the discrepancy between the various reported magic angle analyses. 
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A. General expression for anisotropic EELS 

A general expression for anisotropic EELS can be worked out by substituting the result of Eq. fTU|> back into 
Eq. Q . This gives the cross-section for the partially angular integrated electron energy- loss spectrum for core electron 
excitation in any anisotropic material systems in terms of their macroscopic dielectric function as: 

^(a ,(3o,d) = J^l^[Z x Irn[Tr(e)} + (£„ - U) £cos Xj cos*^'] (25) 

i,j 

As the trace of the dielectric function 'metric tensor' Tr(e) (=^~] e- 7 - 7 ) is invariant with respect to rotational transfor- 

j 

mation, the first part of the expression does not change with specimen orientation. The second part has a pre-factor 
that vanishes at the magic angle conditions. So the cross-section for core electron excitation using MAEELS is 



B. Physical meaning of the magic angle effect 
Another way to write Ea. l|25[l is as follows: 



As before, the factor (£|| — gives the magic angle condition. However, if we rotate the crystal over all possible 
orientations (O), the averaged value of the angular dependent factors can be written as2&: 



(cosxicosxj) = -S i:j (28) 

This means that the second term in Ea. (|27|l drops out when the cross section is averaged over all the orientation, 
even if the magic angle condition is not satisfied, i.e. 

da 8Trm e f Im[Tr(e)} \ 

dE {ao > M = 1^11 + 3 / (29) 

Thus we demonstrated that the spectrum obtained at magic angle condition is equivalent to that obtained by 
orientational averaging. This reason can be explained mathematically as follows. 

In normal orientational averaging, the orientation can refer either to that of the specimen or the orientation of the 
external perturbation, i.e. the orientation of the applied field defined by the momentum transfer vector q in EELS. 
We can distinguish the averaging over the azimuth angle from to 2n , from averaging over the polar angle form to 
7r. In MAEELS, the azimuth angle averaging is achieved by using an axially placed circular detector. The remaining 
orientational averaging of q over the full polar angle range is not possible to achieve in EELS experiments, but the 
equivalent result may be obtained by integrating over a limited range of polar angles because electron scattering is 
skewed towards the small angle. However, it is not always possible to find the appropriate polar angular range, hence 
the magic angle condition is not trivial. 



C. Comparison with other analysis 

Now we can compare our prediction for the magic angle value with other derivations (see table. 1) and to analysis 
the reasons for the diversity of the values predicted. 

Mcnon and YuanS* and Souche et al22 both derived their values by working out the anisotropic spectral response in 
the uniaxial system, but otherwise their conclusions are identical with ours. Our general result is in agreement with 
their prediction for the specific uniaxial system. 
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Most interesting is that of Paxton et al~ who tried to derive a general theory for the magic angle condition. In 
their paper, the momentum transfer vector q is projected in the laboratory frame (X, Y, Z) defined above, so we have 
the quantum mechanical transitional matrix element as: 

K/|q-r)|i)| 2 = \(f\q x X + q Y Y + q z Z\i)\ 2 

= q 2 x (X) 2 + ql{Y) 2 + q 2 (Z) 2 + 2q x q Y Re[(X) (Y)*} + 2q Y q z Re[(Y)(Z)*} + 2q z q x Re[(Z) (X)*] (30) 

where (X) 2 = \{f\X\i)\ 2 , and {X) = \{f\X\i)\, the (Y) 2 , (Z) 2 , (Y) and (Z) have the similar definition. As discussed 
above, using the weighting given in Eq.lj8|l. the angular integrated cross section can be written as: 

^KA,o)« U ■ ((x) 2 + (Y) 2 ) + e r (Z) 2 (31) 

If £j_ = ^|| = ^o,then we have: 

^(c^/^oc i Q -{(X) 2 + (Y) 2 + (Z) 2 ) (32) 

Paxton et al^i derived the magic angle condition by insisting that the isotropic spectra are given by this equation 
without giving detailed explanation. To show that it is orientation independent, we first consider a case where the 
specimen frame coincides with the lab frame, then 

(X) 2 + (Y) 2 + (Z) 2 = (x) 2 + (y) 2 + (z) 2 (33) 

As the latter is proportional to the trace of the imaginary part of the dielectric function, it is invariable to rotation. 
So Paxton et alsii chose the correct magic angle condition, and should arrive at the same value for the magic angle 
as we do. But for reasons we could not understand, they concluded that the MA|| is about 1.369e- We believe that 
it is a trivial mistake. 

Daniels et al.'s derivation 3 ^ has some similarity with our Ea. pUf) . but they have used the substitution in the beam 
direction coordinate (X,Y,Z) as: 

' qx = qxY cos cj) , g4 . 
Qy = qxv sin 4> 

X = r XY cosV 

Y = rxY sin <p y ' 

where the txy was defined as the magnitude of vector rxv-the position vector of the sample electron in the X-Y plane 
where the collection aperture lies, and qxY has a similar definition. However, Daniels et al assumed that <p ^ <fi', so 
their subsequent calculation can not be correct. 

Zhu et al^ii correctly recognized the importance of the rotationally symmetry of the experimental set-up, i.e. that 
the cross term as shown in Eq.JSJ should vanish in integrating over azimuthal angle, so they focused on estimating 
the weighting of the cross-section along the polar angular range as: 

q\\ = q± (36) 

f d 2 a f d 2 a 

where q~i = l / ( ^^ )rif2. The normalization factor in the denominator is just the equivalent expres- 

sion for the isotropic system and we can define it as N. In the small angle approximation, we have: 

*L BdB WTW) (37b) 

In fact q± = k^O where 9 is the so called mean scattering angle defined in Egerton's bookS, and has a value 
9 = 29e{$o — arctan/3 )/ln(/3o + 1), and q\\ = k 9E- i.e. in his definition 9 — 9e and we can be resolve (3q IA w 1.766e 
according to this relation. For comparison, our equivalent integrals defined in Eq.JSJ) and (II 1|) can be rewritten under 
the small angle approximation as: 

t*=r™w+& (38a) 

A) Q2 
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Thus the guess of Zhu et al is incorrect numerically. By a similar argument, Gloter et ai2£ guessed a different 
weighting of transition with q parallel the beam direction and roughly estimated it with (q-k ) 2 /q 2 fco 2 - They also 
used the isotropic angular distribution as normalized factor. Remarkably, their guess is correct for parallel illumination 
case, but their expression can not correctly account for the convergence effect because it does not consider the scattering 
when the incident and scattered electron beams are not in the same plane as the beam optical axis. 

In summary, we have analyzed the reasons for different prediction of magic angle values, and found all the dis- 
crepancies can be properly accounted for. This suggests that there is no fundamental objections to our theoretical 
model. 

We will discuss the discrepancy between the experimental measurement and the theoretical predicted value of 
magic angle in elsewhere. We want to emphasis here that the discrepancy is not due to errors in the theoretical 
analysis, but rather because of the simplicity of the theoretical assumption or the experimental interpretation. We 
can list a number of factors that might modify the prediction in our model, such as non-dipole transition^, coherent 
scattering effect^, channelling effect, relativistic effeci£^l. But our analysis showed that they may affect the exact 
values of the magic angle, but not the conclusion that magic angle effect, if it exists, applies to all anisotropic systems 
and that the spectra collected under magic angle condition represents an orientation averaging. 



IV. APPLICATIONS 



A. Cross sections in anisotropic systems 



We will discuss special cases for the symmetries exemplified by crystalline systems, although in fact the systems 
can be amorphous or nanostructured. This allows us to write out the explicit form of the dielectric function, hence 
the precise form of the cross-section for the partially angular integrated EELS core loss excitation: 

i) Isotropic systems 

This category includes cubic crystals where the non-diagonal elements vanish and the diagonal elements are identical, 
of which the imaginary parts are set to be Ei. So the cross section for EELS is given by: 



da 



k { ^ + «ii -we cos2 ^ } - ^Sk gi + ^ 



(39) 



dE naoh 2 kn 

where we have used the relation cos 2 Xi = 1~- According to the definition of the reduced integral in Eq.JSJ), the 

i 

factor + 2£j_) is in fact proportional to ln(/3 2 + 1) in the case of parallel illumination, 
ii) Uniaxial systems 

There are hexagonal, tetragonal and rhombohedral crystals in this class where the physical response is unique in 
one direction. The dielectric function is characterized by two different element as follows: 

e" i = j = 3 

e 1 - i = j = \or1 (40) 
ijij 

This means that the cross-section for EELS in this class of system is given by: 

{K±(l + cos 2 X3)+e||sin 2 X3]e2 L + (^||COs 2 X3 + e±sin 2 X3)4} (41) 



dE na§h?k§ 
hi) Orthorhombic systems 

In this case, the dielectric function have three independent diagonal elements e-" , so we have: 



da 87rm P 



dE naoh 2 k 2 



+ -£0 cob 2 fc^l ( 42 ) 



iv) Monoclinic and Triclinic systems 

The off-diagonal elements do not vanish in these systems, so we have: 



da 87rm P 



dE naoh 2 k, 



2 



X^K-L% + tell ~ £-0 cos Xi cos XjVi \ (43) 
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B. 'Magic' orientation in the orthogonal systems 

Here we concentrate on the most commonly encountered orthogonal system. The cross section Ea. (l42|l can be 
rearranged as: 

% - ^ {e } - £fe {<* +«™ + «, - «EC * - ^} («) 

Again the first part is rotationally invariant so it represents the isotropic spectrum. The second part contains 
the information about the specimen rotation. The first bracket represents the factor responsible for the magic angle 
condition. The interesting point is the existence of other factors (cos 2 Xj ~ 1/3) inside the summation over j. This 
suggests that the second part will also vanish if all the brackets within the summation sign equal to zero. They 
uniquely define a specific specimen orientation (\j — 54.7°) which we may label as the 'magic' orientation. Again the 
spectrum so obtained equals that obtained at the magic angle condition or by through orientational averaging. This 
is understandable as cosxj IS the projection of the jth basis vector at the optical axis, so each principle symmetry 
electronic excitation contributes equally. By rotation symmetry about the beam direction, the same result holds for 
the more general case involving convergence illumination. 

In uniaxial systems whose dielectric function has only two variables £|| and e± , the angular integrated cross section 
as shown in Ea. l|41|l becomes 

da &Trm e \ . . . ei + 2e* _ , . , 1 , n , | 

dE = + *«> JL T JL + - ^ )(C ° S * 3 - 3 )(£ " ~ ^ j (45) 

Clearly, the 'magic' orientation defined above is reduced to a 'magic angle' between the z-axis of the sample and 
the optical axis, i.e.%3 = 54.7° in uniaxial system. However, we have distinguished this 'magic angle' of specimen 
orientation with the magic angle for the beam convergence and collection in MAEELS. 

In summary, the 'magic' orientation can provide a set up at which the spectra should be the same as the spectra 
gained at the magic angle for the system where the symmetry is higher than orthogonal, and this 'magic' orientation 
will lose its meaning in a system whose dielectric function has the non-zero off-diagonal elements. 

Our analysis suggests that a better way to represent the anisotropic response of EELS is to write it as a linear 
combination of the orientationally averaged (also called 'isotropic') spectrum and an orientation dependent spectrum. 
In uniaxial systems, the orientation-dependent spectrum can be further expressed as a product of the magic angle 
factor, the magic orientation factor and a dichroic spectrum: 

^Anisotropic = (2£_L 1 ^Average + (£|| — £,±){ c OS 2 X3 ~ g) ' £2 \Dichroic (46) 

where 

I ^2 I An \ 

£2\ Average = ^ (,47a) 

Sl\Dichroic = 4 _ £ 2 ( 47b ) 

This formula should offer a practical way to study the anisotropy in the core electron excitation as well as encoding 
the magic angle and magic orientation conditions. 

C. Connections between MAS for EELS, for NMR and for XAS 

The magic orientation effect has a direct analogue with the 'magic angle effect' in XAS experiments and less directly 
with magic angle spinning nuclear magnetic resonance (MAS-NMR) . 

In surface extended x-ray absorption fine structure (SEXAFS)&£, one can study the bonding length as well as local 
coordination number of the excited atoms. However, contribution of more than one shell to the measured EXAFS 
can result in a polarization dependent measured distance from absorbing atom to the neighboring atoms and affect 
the effective coordinate number. If there is higher than twofold symmetry around the surface normal, the correct 
distance and the real coordinate number can be directly measured if the angle, between the electric field vector E 
of the incident x-ray and the surface normal of the single crystal, is equivalent to 54.7° exactly. This is because the 
system being probed is effectively an uniaxial system. 
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Another famous example is the MAS-NMR technique for solidii^iSk. If the material, whether it is a single crystal, 
polycrystal or power, spun with high speed about an axis which is respect to the applied magnetic field with 54.7° , 
the NMR result will be independent of the orientation of the sample. MAS-NMR is also related to the orientational 
magic angle effect because one is effectively using spinning to create an effective 'uniaxial system' out of powered 
samples. 



V. CONCLUSION 



We have presented a general model describing anisotropy of the core-level electron excitation in EELS measurement 
and have determined the magic angle condition at which the sample orientation becomes irrelevant. After comparing 
our derivation with reported theoretical efforts, we can explain all the reasons for disagreement in the literature 
predicting the value of the magic angle and showed that the differences in no way invalid our approach. Furthermore, 
for the first time, we showed that the magic angle condition is applicable in all anisotropic systems and that the 
spectrum at the magic angle condition is equivalent to the rotational average of the sample. The same analysis can 
also give the general expression for electron energy loss spectroscopy of core electron excitation in anisotropic system. 
In high symmetry cases, it leads to the discovery of the magic orientation condition. Its relation with other 'magic 
angle effect' is clarified. In addition, the analysis shows that EELS in uniaxial system can be written as a sum of the 
effective 'isotropic' spectrum and the linear dichroic spectrum. 
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APPENDIX A: THE GEOMETRICAL EXPLANATION FOR EQ.(l22l) 



The perpendicular component of the momentum transfer q can be seen as a sum of the perpendicular components 
of the initial and the final wave-vector in the case of the convergence beam as: 

q ± =k ± -k/ (Al) 

According to the vector relations shown in fig^ under the small angle approximation, we have: 

q± = k 9 (A2) 
fco^ = k a (A3) 
k f x = k (3 (A4) 

Thus the scattering angle 9 can be written following the vector combination rule as: 

6 = a-[3 (A5) 

where the directional properties of these angular vectors are defined to be the same as the perpendicular components 
of their corresponding wave-vectors. According to the law of cosines, the magnitude of 6 therefore can be written as: 
Eq.(|22J. 

6(a, /3,<j>- ip) 2 = a 2 +f3 2 - 2a/3 cos(0 - 92) 
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TABLE I: The predicted values of the magic collection angle for the parallel illumination case. 



Authors 


Zhui^ 


Menorj2£ 


Paxton2i 


SoucheS 


Daniels^ 


MA||(Af A /0 B ) 


1.8 


4 


1.36 


3.97 


1.98 



I The optical axis 
I 




FIG. 1: In an inelastic scattering, the momentum transfer vector q is determined by the initial and final wave-vector ko and 
k/. The angle between the two wave vectors is defined as the scattering angle 9. In the parallel illumination case, 6 always 
equals to f3 defined as the angle between the wave vector of scattering electrons kf and the optical axis. 




(X,Y,Z) (x,y,z) 

i 

FIG. 2: The projection of the momentum transfer vector q in the sample frames (x,y,z) and in the laboratory frame (X,Y,Z). 
The components qi (qi') are equivalent to qe^ (q-ej'), where the ei (e%') are the basis vectors of the reference frame. Two sets 
of components of q are related by the transformation matrix R whose elements are defined in terms of dot product of the basis 
vectors of the two reference frames. 
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FIG. 3: The magic angle is a function of the Be in the parallel illumination case under the small angle approximation, i.e. 
f3™ A =3.976 E (dashed), and for cases no limited by the small angle approximation(solid), i.e. solution of Eq. llH and 1181 . 
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FIG. 4: Schematic diagram of the inelastic scattering under the small angle approximation (a) and its projection in the plane 
perpendicular to the optical axis (b). The convergence angle effect has been taken into consideration. 
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FIG. 5: The contour expression for magic angles condition in electron energy loss spectroscopy of anisotropic core electron 
excitation. 



